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From quark-gluon plasma to hadron spectra 
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Results on initial transverse energy production based on NLO perturbative QCD cal- 
culation with final state saturation of produced minijets are used to fix the initial energy 
density of produced matter. Assuming rapid thermalization, this provides the initial 
^ ■ conditions for a hydrodynamic description of the expansion of final matter. Given a pre- 
scription of the the decoupling of particles from the thermal system to free particles, final 
transverse spectra of hadrons and integrated quantities like multiplicity and transverse 
energy can be calculated in the central rapidity region. Results are reported and compared 
with measurements. 

1. INTRODUCTION 
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Conservative estimates from the measured transverse energy and particle multiplicity at 
RHIC indicate that energy and particle densities well above those of normal nuclear matter 
are formed in nuclear collisions. Results from perturbative QCD (pQCD) calculation of 
minijet production with final state saturation [IJ| follow well the measured multiplicities. 
However, the calculated transverse energy per unit (pseudo)rapidity is larger by a factor 
~ 2.6 than the measured value. For achieving agreement between the calculation and 
the experiment, the evolution of the produced matter must transfer a large fraction of 
the energy from mid-rapidity towards the fragmentation regions. In the hydrodynamic 
description the mechanism for the energy transfer is the work done by the pressure in the 
expansion. Assuming that initially the expansion is mainly in the longitudinal direction 
leads to an asymmetry between the transverse and longitudinal directions and to the 
desired energy transfer. 

The work and the cooling of the matter during expansion will change the momentum 
distributions both in transverse and the longitudinal direction. In this work we assume 
that in the mid-rapidities the longitudinal flow follows the scaling law, v z = z/t (collision 
takes place at t — and z — 0) and that we can ignore the small longitudinal changes of 
densities. With these restricting assumptions we are not able to address the problem of 
longitudinal momentum distributions but we can calculate the final transverse momentum 
distributions when the particles decouple at the end of the thermal stage. Calculation of 
spectra, the folding of flow with thermal motion, is performed using the Cooper and Fry 
prescription. 

The justification for the use of hydrodynamics can be marginal for some parts of the 
produced system but, first, it is clear that the bulk of the matter is very dense leading to 
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numerous secondary collisions and, second, in describing the energy-momentum transfer 
due to these collisions in terms of hydrodynamics the conservation laws are correctly 
satisfied. 

2. INITIAL CONDITIONS FROM pQCD MINUETS 

The perturbative QCD calculation of minijet production is a momentum space calcu- 
lation. In order to define the initial densities a connection between the momentum of 
the minijet and its space-time formation point is needed. At collider energies the hard 
partons of the colliding nuclei are Lorentz contracted to a region of order 2R A / / y cm << 1 
fm. We consider the collision region as a point in the longitudinal direction and assume 
that the rapidity of the minijet coincides with the space-time rapidity of the formation 
point, y = r\ = (1/2) ln[(t + z)/(t — z)\. The formation (proper) time we take to be the 
inverse of the saturation scale, tq = l/p sat - Thus the minijet matter forms along the 

hyperbola t = \J ' z 2 + Tq with initial longitudinal flow velocity v z (tq) = z/t. 

The basic quantities for the minijet production in a nucleon-nucleon collision are 
<7jct(Psat, A) and a jet (E T )(p sat , </s, Ay, A) , the minijet cross section and its first 

moment in transverse energy (momentum) for the rapidity interval Ay, both integrated in 
Pt from px = p S at to infinity 0,0,0 • The number of minijets or the transverse energy in Ay 
in a nucleus-nucleus collision is obtained by multiplying the corresponding cross section 
with the nucleon-nucleon luminosity (with an extra factor 2 for the number of minijets) 
which for a central collision is T AB (0), the overlap function of transverse densities for the 
colliding nuclei: T AB (h) = J d 2 sT A (\h - s\)T B (s) = T AB (b) , T A (s) = J+™dzp A (z,s) = 
T A (s) , where b is the impact parameter and s the transverse coordinate. 

Average densities are obtained by dividing with the volume AV = AzAt = T$Ay ttR\. 
This procedure is easily generalized to local density in the transverse plane of the collision. 
The nucleon-nucleon luminosity for a transverse area element d 2 s is T^(|b — s\)T B (s) and 
the volume element dV = dzd 2 s = rAyd 2 s leading to 

npQCD(ro ' s) = = d^ 2TA(|b " s|)Ts(s)(7jet ' 

and 

dE T 1 

epQCD(T °' S) = ^dWs = ^T A (\b-s\)T B (s)a ]Ct (E T ). 

We next make a bold assumption of fast thermalization: we assume that the thermal- 
ization time scale is the same as that of production and for simplicity we take the system 
to be thermal right at the formation time r . Two simple arguments can be given in favour 
of this assumption. First, if we calculate the temperature either from e(r ) or n(r ), the 
result is very closely the same. This means that for thermalization only collisions causing 
energy and momentum transfer are needed, number changing reactions are not essential. 
Second, the time scale To is associated with the saturation scale and the formation time of 
more energetic (mini)jets could be shorter so that some secondary interactions take place 
already before the overall formation time. 

In Figure 1 the transverse profile of the initial energy distribution is shown for a gold- 
on-gold collision at RHIC (dashed line) and lead-on-lead collision at LHC energy (solid 
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Figure 1. Transverse dependence of 
the initial energy distribution for a 
gold-on-gold collision at RHIC (dashed 
line) and lead-on-lead collision at LHC 
energy (solid line). The average value 
of the saturation scale is p sa t = 1-16 
GeV at RHIC and 2.03 GeV at LHC 
with formation times 0.170 and 0.100 
fm/c, respectively. 



line). We comment on the values of multiplicity and transverse energy of initial minijets 
when discussing later the results on final hadron spectra. 

3. HYDRODYNAMICS OF EXPANSION 

We will assume isentropic expansion. Studies of viscous effects indicate that the entropy 
production during expansion does not produce significant effects 0. The longitudinal 
flow is taken to be boost invariant and to scale at central rapidity. Then v z = z/t 
or equivalently y = rj and the energy density and transverse flow velocity are of form 
e = e(r, r) and vt = v T (r,r) J/]]. With cylindrical symmetry the equations reduce to 
1+1-dimensional form and finding out the solutions numerically is straightforward. 

To close the set of hydrodynamical equations the equation of state (EoS) is needed. We 
have assumed an ideal QGP at high temperatures and a hadron gas including all hadrons 
and hadron resonances up to the mass 2 GeV. Repulsion of hadrons is described by mean 
field with K = 450 MeV and to induce a phase transition a bag constant is included in the 
EoS. Calculations presented here are performed with the bag constant B l l A = 235 MeV 
leading to value T c = 165 MeV for the phase transition temperature f|||. 

Finally, to turn the hydrodynamic quantities, the transverse velocity Vt(t, r) and the lo- 
cal temperature T(r, r) to observable quantities, the thermal motion, characterized by T, 
must be folded with the flow motion vt- This is done using the Cooper and Fry pre- 
scription 0. Final particles are assumed to decouple from the thermal phase at a given 
density ed ec or equivalently in the zero baryon number case at temperature T dec . The 
condition T(r, r) = Td ec defines a hypersurface and the particle spectra are obtained 
as the net particle number flow through this surface. For the boost-invariant cylindrically 
symmetric case the integration can be reduced to one- dimensional integrals: 

+m T I (n-f r v r ^Y)K 1 (n%. 1 ^)dr]. (1) 
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Figure 2. Thermal spectra exhibit trans- 
verse mass scaling when v r = 0. For non- 
zero transverse velocity (v r = 0.6) trans- 
verse spectra of particles depend not only 
on Mt but also on M. Temperatures 
are arbitrary but so chosen that all spec- 
tra have approximately the same slope at 
large Mt- 



Without transverse flow the first term is missing and the second term depends only on 
the transverse mass Mt = \Jp"r + M 2 and temperature T. Equation [I] shows explicitly 
the breaking of Mt scaling through the pt dependent terms when Vt ^ 0. This is shown 
in Figure 2 where pion, kaon and proton spectra are shown for Vt — and Vt = 0.6. 
The temperatures, T = 220 MeV when Vt — and T = 120 MeV when vt = 0.6, are so 
chosen that the slopes at large tut are similar. For v t = the mass shows up only in 
the starting point of the spectrum. (There is a difference between bosons and fermions 
from the alternating sign for the latter in the summation.) For non-zero transverse flow 
the spectra deviate from the scaling behaviour mainly in the region pt<M. When pt is 
clearly bigger than m, the spectra become similar since pt ~ Mt- In nuclear collisions 
the spectra get contributions over a range of transverse velocities, but the general features 
remain the same as in Figure 2 for vt = 0.6. 

The buildup of transverse flow is one manifestation of the collective behaviour of a 
dense system. The estimates, using as input the Lorentz contracted nuclear geometry 
and the observed multiplicities, give quite high particle densities. With such densities 
it would be difficult to understand the dynamics of the produced matter if no collective 
effects are observed since this would indicate that the final state particles interact very 
weakly. 

In any approach based on dominance of hard processes in the initial particle production 
the predicted transverse energy will be consistent with observations only if the system 
behaves collectively: the initially produced transverse energy is much higher than the 
observed one as discussed in detail in Section 4. In the fast longitudinal hydrodynamic 
expansion a large fraction of this energy is transferred into the longitudinal motion which 
should show up in the broadening of rapidity spectra from their initial shape. In the 
present boost invariant calculation only the energy loss in the central rapidity region can 
calculated. 
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4. RESULTS AND COMPARISON WITH EXPERIMENT 

The calculations are performed for central collisions when the produced matter distri- 
butions are cylindrically symmetric. Since the RHIC data is given for a 6 % centrality 
cut, we use in calculating the initial minijet production an effective mass number for nu- 
cleons fixed to equal the average number of participants for that centrality cut; for details 
see ||. We assume a similar centrality cut in predictions for the LHC energies. To show 
the magnitude of the change from imposing the cut, we perform the calculations also with 
mass numbers A = 197 for gold at RHIC and A = 208 for lead at LHC. 
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Figure 3. Charged particle multiplicity 
as function of collision energy p|. Calcu- 
lations (triangles and squares) are for the 
full (open) and effective (filled) mass num- 
ber determined from centrality. 



Figure 4. As in the left panel but for total 
transverse energy ||. The initial trans- 
verse energy at =130 GeV is also dis- 
played. 



I will begin by presenting results on p-r~integrated quantities and finish by discussing 
the transverse spectra. All results are for mid-rapidity since the minijet calculation is 
performed at y — and boost invariance is assumed in the hydrodynamic calculation. 

Figures 3 and 4 display the dependence of the charged particle multiplicity and the total 
transverse energy on the collision energy. The charged particle multiplicity has been mea- 
sured by PHOBOS Collaboration [[□]] at three different collision energies, ^/snn = 56, 130 
and 200 GeV. There are some uncertainties in the normalization of the minijet calcula- 
tion. E.g., the saturation condition for the cut-off momentum in minijet production, 
p sat , could contain extra factors of order 1 but the simplest choice with no such factors 
is taken. However, once the normalization is fixed, the energy dependence is strongly 
constrained. There are still uncertainties common to all hard scattering approaches like 
those coming from higher order contributions. For the minijet calculation used here for 
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Figure 5. Effective temperatures determined as inverse slopes of pion, kaon and proton 
spectra at the LHC energy y/s = 5500 GeV for two different transverse mass intervals. 



the initial conditions, the contribution from the NLO terms has been calculated for the 
transverse energy; for details see [fl , 10f| . The agreement of calculated multiplicities with 



the PHOBOS data is very good both in normalization and the energy dependence. 
The values measured at y/s by PHENIX and STAR |l3j Collaborations are 5 ... 10 
% higher. 

For transverse energy the value measured by PHENIX Collaboration at ^/s = 130 GeV 
[[lljj , 578 GeV, is ~ 6 % smaller than the calculated value 614 GeV. This should be 
compared with the initial value ~ 1550 GeV, also depicted in Figure 4 showing an almost 
a factor of 3 reduction of (IEt/cIt] due to the expansion. On the other hand, as pointed 
out above, the multiplicity measured by PHENIX [12 1 is somewhat above our calculated 
value. This is an indication that, even though the energy in the transverse degrees of 
freedom is reduced by a large factor during the expansion, our calculated spectra, to be 
discussed below, are still slightly shallower than the measured ones. 

Figure 5 shows the results on effective temperatures T e g, the inverse slope of the spec- 
trum, for pions, kaons and protons. When the flow destroys the scaling in transverse 
mass, the slope of the distributions becomes sensitive on the mass of the particle and 
changes rapidly with (or p?) especially in the region My< (2 . . . 3)M as was seen in 
Figure 2. For this reason T e fj is displayed for two different transverse mass intervals, 
0.5 < M T -M < 1.5 and 2.5 < M T -M < 3.5, in the left and the right panels of the 
Figure, respectively. At the smaller transverse mass values the mass dependence of slopes 
is very strong. 

At y/s = 130 GeV, the STAR collaboration has data on inverse slopes at small trans- 
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verse mass, 0.05 GeV < M^ — M < 0.45 GeV | 15fl . These are displayed in Table 1 (in 
MeV units) with inverse slopes from the calculated spectra at M^ — M = 0.3 GeV. One 
observes that the calculated values are slightly larger than the measured ones but the 
ratios of the measured or the calculated slopes for different particles are quite similar. 





Pion 


Kaon 


Proton 


Measured 


190 


300 


565 


Ratios 


1 


1.58 


2.97 


Calculated 


238 


387 


652 


Ratios 


1 


1.63 


2.74 



Table 1 

Inverse slopes, T e g in MeV, of pion, kaon 
and proton spectra measured at RHIC at 
yfs = 130 GeV [15] and compared with 
calculated values. 



Figure 6 shows the calculated average transverse momentum for different particles. In 
addition to the STAR measurement |5| at RHIC, also the NA49 [1§] and the UAl JT?| 

data points are depicted. Comparable to the STAR measurement is (pr)aii = J (Pt) c& 1~ 
culated for all particles and shown as the filled circle. The calculated value is somewhat 
too large as were also the inverse slopes. Comparison of the STAR measurement with the 
UAl point shows clearly a strong nuclear effect and the increase from the NA49 measure- 
ment indicates that the effect grows with energy. In a thermal model with expansion, this 
behaviour is understood in terms of increasing collective flow as the initial densities grow 
with the collision energy. 
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Figure 6. The average pt of pions, kaons and 
nucleons as a function of collision energy. In 
addition to the STAR jfj] data point, also NA49 
[16] and UAl [O results are shown. 
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Transverse spectra of pions, kaons and protons are shown in Figure 7 at y/s = 200 GeV 
and 5500 GeV. The breaking of Mt scaling is clearly seen. It should be emphasized 
that the normalization of kaons and protons relative to pions depends strongly on the 
decoupling temperature which is here assumed to be the same for all particles, Td ec — 120 



MeV. It is more likely that the number of kaons |18| and protons freezes out earlier than 
at the kinetic freeze-out. This would increase their normalization relative to pions but the 
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shape of the spectrum could change less since the elastic collisions could still take place. 
In calculating the strange particle spectra, full chemical equilibrium is used at freeze-out. 
At lower energies strangeness is known to deviate from the equilibrium values We 



have also ignored the non-zero net baryon number which affects the normalization of 
(anti)proton spectrum at RHIC but less at LHC. 




M T [GeV] M x [GeV] 



Figure 7. Transverse mass distributions of pions, kaons and protons. In calculating the 
proton spectrum, zero baryon chemical potential is used. The spectra contain also the 
particles from decays of resonances up to £(1385). 



In Figure 8 the calculated transverse momentum spectrum of negative particles at 
y/s = 130 GeV is compared with the spectrum measured by the STAR Collaboration 
JT3] . There is an indication of difference in shape at small transverse momenta. In the 
calculation the turning down of the spectrum comes from the Jacobian, pt/Mt- At large 
values of transverse momentum the calculated spectrum is somewhat above the measured 
one. This was already seen when comparing the values of inverse slopes and average 
transverse momenta. In terms of the total reduction of energy in the transverse degrees 
of freedom by almost a factor of three this is a small effect and we have not tried to fine 
tune the calculation. 



5. CONCLUSIONS 

We have combined a minijet calculation of initial particle production at collider energies 
with a hydrodynamic treatment for describing the expansion of the produced matter. This 
allows us to predict particle spectra, particle multiplicities and the transverse energy of 
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Figure 8. Transverse momentum spec- 
trum of negative particles. The data 
points, presented as filled circles are from 
STAR Collaboration Jl5[ . The calculated 
spectrum (thick solid line) contains the 
contributions of 7r~, K~ and antiproton, 
shown also separately. 



final particles as a function of the mass number of the colliding nuclei and the collision 
energy. 

The prediction from the minijet calculation of the initial particle production which is 
most robust against details of the treatment of the expansion is the multiplicity density, 
dN/dr] (or dN/dy), since this assumes only an approximate conservation of entropy. Since, 
as was shown in Section 2, the energy and the number density of the initially produced 
minijets correspond to the same thermal state, the number changing reactions are not 
important for thermalization and the entropy production after initial minijet production 
can be expected to be small. The calculations agree very well with the PHOBOS data at 
each of the three RHIC energies. At 130 GeV the PHENIX and STAR results are above 
the calculated value but they are also higher than the PHOBOS result. 

The initial energy per produced parton is quite large. E.g. the calculated average trans- 
verse momentum of initial minijets is (pr)i — 1-47 GeV at the RHIC energy yfs = 130 
GeV and 3.2 GeV at the full LHC energy yfs = 5500 GeV. After the expansion stage 
the calculated value at yfs = 130 GeV is (pr)/ = 0.594 GeV to be compared with the 
measured value 0.514 GeV. This slight discrepancy in the predicted and measured trans- 
verse momentum quantities shows up also in the inverse slopes or effective temperatures 
determined from the spectra. E.g., the calculated values of T e g at Mt = 0.3 GeV are 
20 ... 25 % larger at yfs = 130 GeV than the measured ones. However, the predicted mass 
dependence comes out qualitatively right. The thermal model with collective transverse 
flow predicts a strong dependence of the slope both on mass and transverse mass in the 
transverse mass range Mt — M<2M. 

The results from the minijet calculation with saturation assumption combined with 
the hydrodynamic description of the expansion of final matter look quite promising when 
compared with the first measurements at RHIC. This description of nuclear collisions 
provides a well constrained framework for the study of other signals, like the photon and 
lepton pair emission or the evolution of strangeness and the heavier flavours during the 
final dense state. 
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